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A novel high-throughput second-order-correlation measurement system is developed which records 
and makes use of all the arrival times of photons detected at both start and stop detectors. This 
system is suitable particularly for a light source having a high photon flux and a long coherence time 
since it is more efficient than conventional methods by an amount equal to the product of the count 
rate and the correlation time of the light source. We have used this system in carefully investigating 
the dead time effects of detectors and photon counters on the second-order correlation function in 
the two-detector configuration. For a non- stationary light source, distortion of original signal was 
observed at high photon flux. A systematic way of calibrating the second-order correlation function 
has been devised by introducing a concept of an effective dead time of the entire measurement 
system. 

PACS numbers: 



I. INTRODUCTION 

A second-order correlation function is an intensity- 
intensity correlation function, having information on 
both photon statistics and dynamics of the light genera- 
tion process of a light source. It was first introduced by 
Hanbury Brown and Twiss in order to measure the angu- 
lar separation of binary stars ^ , and later it was applied 
to property measurement of various light sources such as 
measuring the coherence time of thermal light 0, get- 
ting information on the nature of scatterers [3] and sur- 
veying the correlation properties of laser light near laser 
threshold 4]. More recently, it was used in probing the 
nonclassical nature of light such as antibunching 5] and 
sub-Poissonian photon statistics @ . 

Much effort has been made in devising a precise and ef- 
ficient apparatus to measure the second-order correlation 
function. The first successful time-resolved measurement 
was done by using a single detector, a single variable 
delay generator and a coincident circuit to measure the 
coherence time of low-pressure gas discharge in a single 
^^^Hg isotope |2l| . This technique has a limitation in get- 
ting the correct correlation function near zero time delay 
due to the imperfectness of the detector such as spurious 
multiple pulse emission and incapability of detection for 
a finite amount of time just after detecting real photons. 
These effects are referred to as after-pulsing and dead 
time effects, respectively. 

To overcome this limitation, a two-detector configura- 
tion has been adopted in which a light beam is divided 
into two parts and two photodetectors are used to record 
photons arrived at each detector 0, H, 0, . Unlike 
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the single-detector configuration, the spuriously emitted 
photons and dead photons at one detector are completely 
uncorrelated to the photons detected on the other detec- 
tor so that the contribution of after-pulsing and dead 
time effects are equally spread over the entire measure- 
ment time. This allows measurements to be extended 
down to the zero time delay. 

More sophisticated correlators, which made use of mul- 
tiple time delays, were developed which digitized the time 
interval between a start pulse from one photodetector 
and multi ple stop pu lses from the other photodetector 
at a time [ll|, [12, • The number of pulse pairs corre- 
sponding to a given delay is registered on a corresponding 
counter and the result is proportional to the second-order 
correlation function. Such a device is called a multi-stop 
time-to-digital converter (MSTDC). This method is more 
efficient than previous method owing to the multiple de- 
lay generators. 

Even though a two-detector configuration can effec- 
tively remove the artifacts caused by the detector imper- 
fectness, this benefit is true only for a stationary source, 
the intensity of which is independent of time. Since the 
probability of spurious emissions or losing photons due 
to the detector dead time depends on the intensity, the 
measured intensity profile can be greatly distorted for 
a non-stationary source. In other words, the spurious 
emissions or lost photons are not equally spread over the 
entire measurement time but dependent on time. This 
can cause distortion in the second-order correlation mea- 
surement. 

In the present study, we have carefully investigated 
the limitation of the two-detector configuration for non- 
stationary sources and have devised a systematic way 
of calibrating the second-order correlation function for 
the first time to our knowledge. For this study, we have 
developed a novel second-order correlation measurement 
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system which records all the arrival times of photons de- 
tected at start and stop detectors and make use of all 
the photons detected at the start detector as triggers. 
There are no waste of photons at the start detector and 
thus our system can be much more efficient than MSTDC 
for a light source having a high photon flux and a long 
coherence time. 



II. EXPERIMENTAL APPARATUS 

A schematic of our second-order correlation measure- 
ment system is shown in Fig.^ Photons are detected by 
two avalanche photodiodes (APDs, model SPCM-AQR- 
13 by PerkinElmer). One detector (APDl) serves as a 
start detector while the other (APD2) as a stop detec- 
tor. Each APD has a dark count rate less than 150 cps 
and a dead time of about 50 ns and generates electrical 
pulses whenever the photons are detected with detection 
efficiency of about 50%. APDs are electrically connected 
to the counter/timing boards (Counter 1, 2) installed in 
two computers (Computer 1, 2). 

Relatively low-cost and commercially available 
counter/timing boards (PCI-6602 by National Instru- 
ments) are used to record the arrival times of the 
electrical pulses and to store them to the computers. 
Each board has its own internal clock of 80 MHz so that 
the time resolution is 12.5ns. Two independent boards 
installed in separate computers are used to prevent the 
crosstalks between boards or computers. 

The internal clock of each counter/timing board has a 
limited accuracy and also has a drift of 50 ppm as the sur- 
rounding temperature changes. The accurate frequency 
of the internal clock in each board has been calibrated by 
counting the arrival times of reference pulses from func- 
tion generator(DS345, Stanford Research System). The 
frequency difference between boards is typically several 
tens of ppm and is accounted for in getting absolute ar- 
rival times from the measured ones. 
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FIG. 1: Schematic of experimental setup. ECDL:Extended 
cavity diode laser(2010M:Newport), A;Acousto-optic modula- 
tor(Isomat 1206C), FG:Function generator(DS345, Stanford 
Research System), APD1,2: Avalanche photodiode(SPCM- 
AQR-13, PerkinElmer), Counter 1,2: NI 6602 counter/timing 
board. Counter/timing boards are installed in computer 1 
and computer 2 separately, and are simultaneously armed by 
a trigger signal generated from computer 0. 



To make the two counters start to count at the same 
time, an additional control computer (computer 0) is 
used to generate a trigger signal to simultaneously arm 
the counters. It has a board with analog outputs and 
digital inputs/outputs (NI6703, National Instruments), 
which can send a TTL signal as a trigger. 

All the arrival times of photons detected on both de- 
tectors are the relative times with respect to the same 
origin defined by the trigger. All the detected photons 
on one detector therefore can be used as start pulses 
with respect to those of the other detector. For this rea- 
son, we call our apparatus a multi-start multi-stop time- 
to-digital converter (MMTDC) compared to the conven- 
tional multi-stop time-to- digital converter (MSTDC). 

MMTDC makes use of all the photons detected at a 
start detector whereas the MSTDC makes use of a single 
photon detected at a start detector as a trigger and mea- 
sures the relative arrival times for a time interval which 
is several times (let us say n times) longer than a correla- 
tion time Tc of a light source. It starts over and repeats 
the next measurement cycle using another single photon 
detected at the start detector after nTc- If the incoming 
photon flux to the start detector is 7, only one photon 
out of n'yTc photons is used in the measurement. There- 
fore, our MMTDC has an efficiency n^Tc times higher 
than that of MSTDC. 

MMTDC is specially advantageous for a light source 
which has a high photon flux and long correlation time 
but has a limited operation time with its second-order 
correlation signal embedded in a large background. The 
microlaser was a good example to fit this cate- 

gory. It had an output photon flux of about 3 Mcps and a 
correlation time of about 10 fj,s such that our new method 
was about 30 times more efficient than that of MSTDC. 
Because of limited oven life time, full time measurement 
could give a signal-to-noise ratio of about 3 even when 
MMTDC was used. We could have obtained a signal-to- 
noise ratio of only 0.55 if we had used the conventional 
MSTDC method, only to fail to resolve the signal. 

Due to limitation in computer memory, the number 
of arrival times recordable at a time is limited by about 
one million counts in our MMTDC setup. To get an 
enough number of data, measurements should be done in 
a sequential way. All computers are connected by eth- 
ernet connections in such a way that they can send and 
receive messages among themselves. The counting com- 
puters (computer 1 and 2) send a message to the control 
computer to notify the end of counting whenever they 
complete a specified number of counting and recording. 
After checking that the control computer has received 
the message, both counting computers prepare a next 
measurement. When the control computer receives the 
message, it sends triggers to the counting computers to 
initiate counting again. The number of sequences is de- 
termined so as to get an enough signal-to-noise ratio. 

To get the second-order correlation function, a his- 
togram is constructed for the time differences between 
all possible pairs made of one of the photon arrival times 
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at start APD and another one at stop APD. To save cal- 
culation time, only the pairs of photons the time differ- 
ence of which are within a certain time window, typically 
chosen 10 times larger than the correlation time of the 
source, arc included in the calculation. The second-order 
correlation function can then be obtained from the nor- 
malization of this histogram by the averaged histogram 
value for much longer delay times than the correlation 
time. 



III. EXPERIMENTS WITH STATIONARY 
SOURCES 



We have tested our MMTDC system using an extended 
cavity diode laser (2010M, Newport) operating far above 
threshold as a test source. The photon statistics of its 
output is supposedly Poissonian and the second-order 
correlation function is thus unity for all delay times. Pho- 
ton flux was 3 Mcps for each detector and the number of 
sequences was 1,000 with 500 kilo counts per sequence. 

A correlation function obtained from the photons mea- 
sured at a single detector exhibits the effect of detector 
dead time and after-pulsing. A figure |2Ia) shows a mea- 
sured result. Two types of dip below 1 appear and the 
one near the zero time delay corresponds to the detector 
dead time. This can be confirmed from the output pulse 
shape of APD shown in Fig.|2Ic). The full width of the 
pulse was measured to be about 50 ns. The stop photons 
following a given start photon within this time window 
are completely ignored and thus a dip with a depth of 
unity appears. 

Another dip is extended to 250 ns and its depth is 
about 0.25, which means that the following stop photons 
are partially ignored with a probability of 25%. This 
dip originates from the counter/timing board. Usually 
counter electronics also have dead times since it takes a 
finite amount of time to record measured arrival times. 
In our case, since there is no onboard memory in the 
counter/timing board, it has to transmit the arrival times 
to the computer memory through DMA (Direct Memory 
Access). Since the data transfer rate through DMA is 
limited by 100 MBps, some counts can be missed if the 
time interval between photons is too short to be trans- 
ferred. The probability to miss a count depends on the 
time interval between successive photons. We call this 
loss of counts the dead time effect of the counter as an 
analogy to the detector dead time. 

The correlation function obtained from the photons 
detected at two detectors is shown in Fig. Efb). The 
effect of the detector and counter dead time is completely 
removed and the value of the correlation was very close to 
unity for all delay times, as expected. Since the number 
of photon pairs per bin amounted to 10®, the standard 
deviation from unity was only lO"''. Utilizing all the 
photons detected at start detector helped to reduce the 
background noise greatly. 
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FIG. 2: (a) Second-order correlation function for a Poissonian 
light source measured in the single-detector configuration, (b) 
The same measured in the two-detector configuration, (c) 
Shape of a single-photon pulse from APD. 



IV. EXPERIMENTS WITH NON-STATIONARY 
SOURCES 



Using our MMTDC second-order correlation measure- 
ment system, we have measured the second-order correla- 
tion function for a non-stationary light source. The out- 
put beam from an extended cavity diode laser was mod- 
ulated by an acousto-optic modulator. The amplitude of 
a driving RF (radio frequency) field to the acousto-optic 
modulator was sinusoidally modulated using a function 
generator so that the intensity of the first-order diffracted 
beam was sinusoidally modulated. Its functional form 
can be written as I(t) = a siii{ujt + (p) + b, where a and b 
have the units of count per second for the photon count- 
ing measurement. Since the response time of AOM is 
measured to be 130 ns, the modulation frequency can be 
safely set up to 1 MHz. For the present experiments, the 
modulation frequency was set as 100 kHz. 
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The normalized second-order correlation function for 
a classical source, applied to the intensity of electromag- 
netic field I{t), is given by 



/2)(r) = 



(1) 



where (...) denotes a time average. The second- 
order correlation function for the sinusoidally modulated 
source is thus calculated to be 



g(^)(r) = 1 + — coswT. 



(2) 



Figure |3f a) and (b) show measured second-order correla- 
tion functions. Mean count rates for individual detectors 
were about 0.6 Mcps and 1000 sequences of measure- 
ments were done with each sequence counting 300 kilo 
counts for each detector. It took about 1000 seconds in- 
cluding sequencing procedure. The adjacent points were 
added up and thus the time resolution was 125 ns. 

Figure|3Ja) is the second-order correlation function ob- 
tained from photons detected at a single detector and 
thus shows a sharp dip near time delay zero. The dip re- 
sults from the dead times of the detector and the counter. 
On the other hand, Fig. E^b) is obtained from the pho- 
tons detected on two detectors, APDl and APD2. The 
central dip is completely eliminated, as in the case of sta- 
tionary light sources. The normalized shot noise is only 
about 0.06% due to the large mean counts per bin of 
about 2.8 million. The contrast ratios a^/26^ of single- 
and two-detector configurations are almost the same. 



A. Dead Time Effect on Non-Stationary Intensity 

Even though the two-detector configuration can elim- 
inate the detector dead time effect near the zero time 
delay, the detector dead time still can affect the correla- 
tion measurements for non-stationary sources when the 
detector counting rate becomes comparable to the in- 
verse of the detection dead time. Since the probability 
to miss photons is dependent on the intensity, the mea- 
sured time-varying intensity profile can be distorted. 

We used photon counting method to measure the time- 
dependent intensity. Since our MMTDC was based on 
the photon counting method, we used photon counting 
technique rather than photo current measurement for 
a consistent quantitative relation between the intensity 
profile and the second-order correlation function. In or- 
der to obtain enough number of counts per counting bin 
for intensity measurement, which has to be done in a 
fashion of single shot measurement, we used a slow inten- 
sity modulation frequency of 0.1 Hz for a given incoming 
photon flux of around 1 MHz. 

Photons were counted for every counting bin of 0.01 
sec and the counted numbers were transferred at the end 
of each bin. Since the dead time of the counter occurs 
when the time interval between successive data transfers 



is shorter than 250ns, this measurement should be free 
from the counter dead time effect. At 1 MHz counting 
rate, the mean counts N per bin was 10* counts and thus 
the normalized noise \/N/N was 0.01, which means we 
could resolve an intensity modulation whose contrast was 
as small as 0.01. 

We repeated the measurement for various mean count 
rates b while keeping the contrast of the intensity a/6 con- 
stant. To do so, we fixed the radio frequency field driving 
the acousto-optic modulator and varied the mean inten- 
sity using neutral density filters. The intensities mea- 
sured in this way were well fitted by sinusoidal functions 
as shown in Fig. QJa). The contrast ratio in Fig. Qfa) 
is about 0.422, which corresponds to the contrast ratio 
at a mean intensity of 3 Mcps in Fig. EJb), where the 
measured results denoted by circular dots show a linear 
decrease as the mean count rate increases. 

This linear decrease can be explained in terms of the 
correction factor a of the detector which needs to be 
multiplied by the measured intensity to get an original 
intensity. 



1 + T,L 



(3) 



where / and /,„ are the original and measured intensities, 
respectively, in the unit of cps and Td is the dead time of 
the detector. If there exists additional dead times such 
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FIG. 3: Typical result of the second-order correlation func- 
tion for a sinusoidally modulated light source measured in (a) 
single-detector configuration and (b) two-detector configura- 
tion. 
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as the dead time of the counting board, Td should be 
replaced with an effective total dead time including all 
dead times, and this subject will be discussed in details 
in the next section. For T^a ^ 1 and Tdh <C 1, 

/„i(t) ~ {I -TdIit))I{t)~b-Tdb^ + {a- 2abTd) sin ujt 
~ 6(1 - Tdb) [1 + (a/6)(l - Tdb) simut] , (4) 

which shows that the contrast ratio is modified from a/b 
to {l~Tdb){a/b). 

Note that the unmodified contrast a/b, fixed in the ex- 
periment, can be determined from the y intercept of a 
linear fit to the measured contrast ratios and the incli- 
nation of the linear fit corresponds to the dead time Td- 
In Fig^Jb), we get a/b ^ 0.5 and Td — 56 ns, which is 
about 10% larger than our initial estimate based on the 
single-photon pulse shape in Fig.[5fb). This discrepancy 
is due to the finite bin size of 12.5 ns, which results in an 
additional broadening of about 6 ns, a half of the bin size, 
in the effect of the detector dead time in the correlation 
function. 
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B. Dead Time Effect on 5^2' (r) 

Since the intensity profile is distorted by the detec- 
tor dead time, the second-order correlation function will 
be also affected. For a measured intensity contrast of 
a'/6', the expected contrast of the second-order correla- 
tion function is They are denoted by circular 
dots in Fig. O For an unmodified contrast a/6, the de- 
pendence of the measured contrast a'/6' on the mean in- 
tensity 6 is given by (1 — Tdb)a/b. The contrast of g^^-* (t) 
is thus expected to be 



a'V26'- 



[(1 - Tdb) (a/6)]' ~ (1 - 2T,6)aV262 (5) 



The dashed line in Fig.Oindicates the expected contrast 
with Td — 56 ns. 

Using MMTDC, we have measured the second-order 
correlation function under the condition identical to the 
one under which we had measured the intensity except 
for the intensity modulation frequency, which is now set 
at lOOkHz, since the correlation measurement need not 
be done in a single-shot fashion as in the intensity mea- 
surement. The square dots in Fig. [S] arc the contrasts 
obtained by fitting the measured 5*-'-'(t) with a sinu- 
soidal function. It decreases linearly as the photon flux 
increases, but the decreasing rate of the contrast ratio is 
1.54 times larger than that expected from the intensity 
measurement. 

This discrepancy is due to the dead time of the 
counter /timing board, which was absent in the preceding 
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FIG. 4: (a) Intensity profile of a sinusoidally modulated light 
source measured by the photon counting method, (b) Con- 
trasts of the intensity. Circular dots: measured contrast ratio 
from fitting the intensity, line: theoretical calculation includ- 
ing the effect of detector dead time of 56 ns. 



FIG. 5: Contrast of the second-order correlation function. 
Square dots: calculated from the measured g^^Hr) (100 kHz), 
circular dots: expected from the measured contrasts of the 
intensity (0.1 Hz modulation). A solid and dashed lines are 
theoretical results with the effective dead time of 80 ns and 
56 ns, respectively. 
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intensity measurement performed at 0.1 Hz modulation 
frequency with a data transfer rate of lOOHz. In the 
second-order correlation measurement, regardless of the 
modulation frequency, since all of the arrival times are 
recorded, the data transfer rate can be much faster than 
the inverse of 250 ns, which is the maximum dead time 
of the counter. When the time difference of two succes- 
sive photons is shorter than 250 ns, there exists a finite 
chance that the following photon is ignored. A complex- 
ity arises since this chance is not alway unity. It can be 
anywhere between and 1. We thus need to find an effec- 
tive dead time which can properly includes both detector 
and counter dead time effects. 



V. ANALYSIS AND CALIBRATION METHOD 

To understand the dead time effect of counter, we nu- 
merically simulated the effect of a partial dead time. 
A Poissonian light source was simulated using a ran- 
dom number generation algorithm. In the simulation, 
if the time difference between two successive photons are 
shorter than 50 ns, the following photon is omitted with 
a probability Pl of 50%. This probability would be unity 
for the case of the detector dead time. Figure Efa) is the 
second-order correlation function of this simulated source 
and shows a dip with a half width the same as the dead 
time and a depth equal to Pl, 50%. The correction fac- 
tor a is calculated as a function of the mean intensity 
and shown in Fig. |HIb). The result can be well fitted by 
the following relation. 



= 1 + PLTdl, 



(6) 



Note that the linearity coefficient in a with respect to 
the intensity is not Td as in Eq. ((SJ but PhTd- We can 
generalize this observation and expect that the detector 
dead time Td in Eq. (O will be replaced by PhTd- This 
expectation has been confirmed by our numerical simu- 
lations. Therefore, we call Tes = PhTd an effective dead 
time and it is equal to the area of the dip (r > only) 
around the origin in the second-order correlation function 
obtained in the single-detector configuration (Fig. IHIa)). 
In general, we can experimentally determine the effec- 
tive dead time T^^"^^ of an entire detection system from 
the dead time distribution, i.e., the shape of the second- 
order correlation function around the origin, calculated 
with photons detected on a single detector only. 

The dead time distribution of our MMTDC system is 
shown in Fig. Efa) . The total area of the dip near the 
origin is 80 ns, which is our effective dead time 7^,^°*'' for 
the entire measurement system. The area of the central 
dip only equals Td=55 ns, originating from the detec- 
tor dead time. The ratio T^^^^ /Td is 1.45, and thus we 



expect that the inclination of the observed contrast in 
g'^'^^r) with respect to the counting rate b would be 1.45 
times greater than that expected from the intensity mea- 
surement free from the counter dead time effect. This is 
1.2- 
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FIG. 6: Numerical simulation result for a partial dead time 
effect. A dead time of 50 ns and a probability to lose photon 
of 0.5 are assumed, (a) A second-order correlation function, 
(b) Correction factors as a function of the mean count rate of 
incoming photons. Solid line is a linear fit. 



what we have observed in Fig|Sl where the observed ratio 
is 1.54, showing only 6% discrepancy from our expecta- 
tion. 

From these observations, we can establish a calibra- 
tion method for the second-order correlation function of a 
non-stationary light source. We first measure the second- 
order correlation function using the single detector con- 
figuration and determine the effective dead time T^g°*'' of 
the entire system from the dead time distribution around 
the origin. The second-order correlation function is then 
measured again, but this time in the two detector config- 
uration. Since the correction factor in Eq. (6) is indepen- 
dent of the modulation frequency, we can calibrate the 
measured correlation function by deviding its amplitude 
by (1 - 2rjg°*^6), where b is the mean count rate. 

This work was supported by the Ministry of Science 
and Technology of Korea and NSF grant no. 9876974- 
PHY. 
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